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Abstract The spins often stellar black holes have been measured using the continuum- 
fitting method. These black holes are located in two distinct classes of X-ray binary 
systems, one that is persistently X-ray bright and another that is transient. Both the 
persistent and transient black holes remain for long periods in a state where their 
spectra are dominated by a thermal accretion disk component. The spin of a black 
hole of known mass and distance can be measured by fitting this thermal continuum 
spectrum to the thin-disk model of Novikov and Thorne; the key fit parameter is the 
radius of the inner edge of the black hole's accretion disk. Strong observational and 
theoretical evidence links the inner-disk radius to the radius of the innermost stable cir- 
cular orbit, which is trivially related to the dimensionless spin parameter a* < !)• 
The ten spins that have so far been measured by this continuum-fitting method range 
widely from a* ft; to a* > 0.95. The robustness of the method is demonstrated 
by the dozens or hundreds of independent and consistent measurements of spin that 
have been obtained for several black holes, and through careful consideration of many 
sources of systematic error. Among the results discussed is a sharp dichotomy between 
the transient and persistent black holes; the latter have higher spins and larger masses. 
Also discussed is recently discovered evidence in the transient sources for a correlation 
between the power of ballistic jets and black hole spin. 

Keywords black hole physics • accretion disks ■ X-ray binaries ■ stars: winds, outflows 



1 Introduction 



In his Ryerson Lecture. IChandrasekharl ( 19751 ) described the Kerr solution as the "most 



shattering experience" of his entire scientific life. He found himself "shuddering be- 
fore the beautiful, the incredible fact" that each of the many trillions of black holes 
in the universe is completely described by a single pair of numbers that specify the 
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black hole's mass and its spin Q. In Chandrasekhar's time, 1910- 1995, the masses of 



CyR X-1 and three other stellar black holes had been estimated ijW ebster & Mmdin 



1972 




1992 





19721: ICowlev et al.lll983l : iMcClintock fc Remillardlll986 : Rcmillard et al. 



Bolton|[__ - ' M l 

Shahbaz et al. I ll994 ). Today, accurate dynamica l mass measurements have been 



achieved for more than a dozen stellar black holes (jMcClintock fc RemillardI I2OO6I : 



Ozel et al.ll2010l : IOrosz et al.ll2011al: ISteeehs et allboiiT as well as for s everal super 



massive black holes, e.g., Sgr A* dChez et al."'2008': 'Cillessen et al.l|2009|), NGC 4258 



(|Herrnstein et ah. 2005^ . and others i Giiltekin ct al . 2009 , and references therein). 



In 1989, the first p ractical approach to measuring black hole spin was suggested 
bv lFabian et al.l 1 19891 ), namely, modeling the relativistically-broadened Fe K emission 



line emitted from the inner accretion disk. The first compelling observation of suc h 
a line was reported just two months before Chandrasekhar died ( Tanaka et al.lE995l ). 



Presently, the spins of more than a dozen black holes have been estimated by modeling 
the "reflected" spectrum of an accretion disk, which includes as its most prominent 
feature the Fe K l ine. For a review of this method of measuring black hole spin, we 
refer the reader to iRevnoldsl \201± . 

It was not until 1997 that a new approach to measuring black hole spin - the 
continuum-fi t ting m ethod - which is the subject of this chapter, was pioneered by 
IZhang et al ] (|l997l ). In brief, in applying this method one fits the thermal contin- 
uum spectrum of a black hole's accretion disk to the relativistic thin-disk model of 
iNovikov fc Thornd 11973*) and thereby determines the radius of the inner edge of the 
disk. One then identifies this radius with the radius of the inn ermost stable circula r 
orbit (-Risco)i which is simply related to the spin parameter a* (|Bardeen et al.lll972l '). 



The method is simple: It is strictly analogous to measuring the radius of a star whose 
fiux, temperature and distance are known. By this analogy, it is clear that it is essential 
to know the luminosity of the accretion disk; hence, one must have estimates of the 
source distance D as well as the disk inclination i. Additionally, one must know M in 
order to scale -Risco ^-nd thereby determine a*. 

In 2006, the co ntinuum-fitting method was employed to e stimate the spins of three 



stellar black holes ( Shafee et al.l200^ : lMcClintock et al.ll2006l ). Presently, ten spins have 



been measured using this method (Section|6}. Not only is the continuum-fitting method 
simple, it is also demonstrably robust. For example, there is strong observational and 
theoretical evidence (discussed in Section |4]) that the disk is truncated quite sharply at 
i?igco- Furthermore, there is an abundance of suitable X-ray spectral data for many 
black holes; consequently, for a given black hole one can typically obtain tens or even 
hundreds of independent measurements of spin that agree to within a few percent 
fSection l4.1|) . The one open question for this method is whether the black hole's spin 
is aligned with the orbital angular momentum vector of the inner disk (Section 15. 4|) . 
Meanwhile, a limitation of the continuu m-fitting method is that it is only readily 
applicable to stellar black holes (but see ljoUev et al.|[2009l : [Czernv et "al]|201l|). while 



the F e K method is applicable to both stellar and supermassive black holes ( Revnoldi 
20li ). 

In order to obtain secure measurements of spin using the continuum-fitting method, 
and to establish the reliability of this method, substantial and comparable effort is re- 
quired on three fronts: (1) The selection and fitting of X-ray spectral data to the 
Novikov-Thorne model (in conjunction with ancillary models); (2) testing and explor- 



^ Spin is commonly expressed in terms of the dimensionless parameter a, = cJ/GMc^, 
where J and At are respectively the angular momentum and mass of the black hole. 
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ing extensions of the Novikov-Thorne model via general relativistic magnetohydrody- 
namic (GRMHD) simulations; and (3) obtaining accurate estimates of D, i and M. 
The first two topics are discussed in Sections [HHS] Concerning the third topic, we refer 
the reader to recent p apers on the measurements of these crucia l parameters for M3 3 
X-7 (jOrosz et al.ll2007ll: LMC X-1 JOrosz et al.ll2009l); A0620-00 JCantrell et al ll2O10|l : 



XTE J1550 -564 JOrosz et alJlboilbl ): Cvg X-1 jReid et al.ll201ll: lOrosz et al.ll2oTTal i 
H1743-322 (|Steiner et al.ll2012al ): and GRS 1915+105 jSteeghs et al.ll2013l ). The 



un- 



certainties in D, i and M are critically important because they dominate the error 
budget in the final determination of a» , including the error incurred by reliance on the 
Novikov-Thorne disk model (Section [5|. 

Initial efforts are under way to use the available continuum-fitting spin data to 
i nvestigate the for mation and evolution of black holes, as well as their host systems 
(IWong et a l."2012^. and to understand how a black hole interacts with its environment 



(e.g., Wang ct al. 200^ : ICooke c t al. 2008). The most important application to date of 
spin data is the discovery of a long-predicted correlation between b allistic jet power and 
black hole spin, which is the subject of Section [T] Very recently, iRussell et al.l (|2013l ) 
challenged the validity of this correlation; Section TTM answers this challenge. 



2 Stellar Black Holes in X-ray Binaries 

There are 24 confirmed black hole binaries : the 23 listed in Table 1 in lOzel et al] 
( 20101 1 plus H1743-322 (jSteiner et all boOQal 'PI. A schematic sketch to scale of 21 of 



these confirmed black-hole systems is shown in Figure [T] 

For decades, it has been customary to define two classes of X-ray binaries, com- 
monly referred to as LMXBs (low-mass X-ray binaries) and HMXBs (high-mass X-ray 
binaries), based o n whether the mass of the secondary star is relatively low or rel- 
atively high (e.g.. iBradt fc McClintoc^ Il983h . Here, we use a different classification 
scheme that differentiates two distinct classes of black hole binaries by the primary 
mode of mass tra nsfer to the black h ole and the effect that this has on the stability of 
the X-ray source (|White et aLlllOOSl l. 



The black holes in five of the 24 systems are steadily fed by the winds of mas- 
sive O-supergiant or Wolf-Rayet companions, and consequently the ir bolometric X-ray 
luminosities are relative ly stable (e.g., see Section 4.2.2.1 in McCl intock fc RemillardI 
I2OO6I : lOrosz et al.ll2007h . Sketches of three of these systems appear in the top-right 
corner of Figure [1] We refer to these systems and their black holes as "persistent." 

The black holes in the remaining systems are fed by Roche lobe overflow through 
the LI point, and all of them have been observed to vary in luminosity by factors of 
> 100 (~ 10* in extreme cases; e.g., see iNaravan fc McClintock l200a ). We refer to 
these systems and their black holes as "transient." 



^ Apart from H1743-322, our selection is based on firm dynamical evidence, and we therefore 
exclude some important systems for w hich there is significant evidence that the primary is a 
black hole, e.g., Cyg X -3 llZdziarski et al. 2013), or a strongpresumption that it is, e.g., SS433 
dBegelman et al.||2006l'l and 4U 1957-1-11 llNowak et al.ll20ia) 
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Fig. 1 Schematic sketch to 
scale of 21 black hole bi- 
naries (see scale and leg- 
end in the upper-left corner). 
The tidally-distorted shapes 
of the companion stars are 
accurately rendered in Roche 
geometry. The black holes 
are located in the center of 
the disks. A disk's tilt in- 
dicates the inclination an- 
gle i of the binary, where 
i = corresponds to a 
system that is viewed face- 
on; e.g., i = 21° for 4U 
1543-47 (bottom right) and 
i = 75° for M33 X-7 (top 
right). The size of a system is 
largely set by the orbital pe- 
riod, which ranges from 33.9 
days for the giant system 
GRS 19154-105 to 0.2 days for 
tiny XTE J1118-I-480. Three 
well-studied persistent sys- 
tems with their supergiant 
secondaries are located in the 
upper-right corner. (Figure 
courtesy of J. Orosz.) 



2.1 Persistent Black Hole Binaries 

These systems are distinguished by the large masses of their secondary stars (20 Mq — 
70 Mq) and by the extreme optical/UV luminosities of these stars, which exceed the 
X-ray luminosities of their black hole companions. Consequently, the effects of X-ray 
heating are minimal and the optical star dominates the optical properties of the system. 
The key distinguishing feature of these systems is their X-ray persistence, which is a 
consequence of the star's massive stellar wind (~ 10^'^ — W^^Mq yr~^), a significant 
fraction of which is captured by the black hole. 

Because the secondaries are massive these systems are obviously young (< lO^yr). 
They are also very rare: There is only one confirmed system in the Galaxy, Cyg X-1, 
and, despite many deep Chandra and XMM-Newton X-ray observations of Local Group 
galaxies, only four other such systems have been discovered, one each in the LMC, M33, 
IC 10 and NGC 300. 

In this review, we do not consider further the two persistent systems that contain 
Wolf-Rayet secondaries, namely IC 10 X-1 and NGC 300 X-1, because the masses of 
their black holes depend strongly on the very uncertain masses of their secondaries, 
and also because no attempt has so far been made to estimate their spins. By contrast, 
the three remaining persistent systems - M33 X-7, LMC X-1 and Cyg X-1 - have 
well-determined values of both mass and spin (see Section [B]). Relative to the black 
hole primaries in the transient systems (apart from GRS 1915+105), the black holes 
in these three persistent systems have large masses, M = 11 — 16 Mq, and high spins 
that range from a* = 0.84 to a* > 0.95, a point that we return to in Section [6. II 
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Fig. 2 Outburst cycle of 
XTE J1859+226 in 1999. 
The dashed line (top three 
panels) marks the time of 
peak radio flux (panel d). 
The 1-day radio spike 
(panel d) is shown fully 
resolved in Fig ure 2 i n 
iBrocksoDP et al.l ||2002| ). 
The red crosses (panel 

b) indicate times when 
the X-ray spectrum is 
dominated by the thermal 
component. These BATSE 
and RXTE/ASM X-ray, 
and Merlin (and other) 
radio data (panels a, 
b and d, respectively) 

appear in Figu re 1 in 

iBrocksopp et al.l ||2002|) 
and the optical data (panel 

c) ap pear in Figure 2 i n 
Sanch ez-Fernandez et al.l 
II200HV For further details, 
consult the references. 



2.2 Transient Black Hole Binaries 

With few exceptions, the 18 transient black hole binaries (hereafter simply referred to 
as transients) manifest and then rise to maximum X-ray luminosity on a timescale of 
several days, thereafter returning to a quiescent state over a period of many tens or 
hundreds of days, as illustrated in Figure [21 The masses of the black holes in these 
systems are relatively low, as are their spins (with the exception of GRS 1915-1-105), 
and their orbital periods range widely from 0.2-33.9 days. By comparison, the orbital 
periods of the persistent sy stems span a relatively narrow range. Th e transients are. 



periods ot the persistent sy stems span a relatively narrow range, ih e 
on average, likely Gyrs old l|White fc Ghoshlll998l : iFragos et al.ll2013l '). 



During a ma jor outb urst, the peak luminosities of transient sources approach the 
Eddington limit (jSteiner et al. 2013), while in quiescence their luminosities are typically 
in the range lO"^'^ to 10"*^ of Eddington paravan fc McClintockl l20(i«i ). Figure i] 
shows X-ray, optical and radio light curves of a typical short-period transient. The 
optical emission is generated largely by reprocessing of X-rays in the accretion disk, 
and the radio outburst is primarily the result of synchrotron emission produced in a 
jet. The ballistic jets, which are the subject of Section [T] are launched very near the 
time of peak radio emission (panel d), which in the case of XTE J1859+226 occurred 
just 0.5 days afte r the X-ray l uminosity (panel b) peaked as the source approached its 
Eddington limit IjSteiner et al.; 2013). Spin can be reliably measured when the thermal 
component dominates the spectrum during the latter part of the outburst cycle (panel 
b). For a complete and state-coded version of the X-ray light curve of XTE J1859-f 226, 
see Figure 8b in Rcmillard fc McClintock (2006). 

There are several oddballs among the transients: Four of them have relatively mas- 
sive secondaries, ~ 2 — 6 Mq, compared to the typical value of < 1 Mq ( Charles fc Cod 
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bood l. GRS 1915+105 has remained very luminous continuo usly since its appearance 
in 19 92, and GX 339-4 never reaches a deep quiescent state (jMcClintock &: RemillardI 
LMC X-3 is almost alw ays active and highly variable fSection l4.ip . although it 
does have extended low states ( Smale &: Bovdir2012h . 



3 The Continuum-Fitting Method 

The two foundations of the continuum-fitting method are (1) the existence of an ISCO 
for a test particle orbiting a black hole and (2) the strong observational and theoretical 
evidence that - for a wide range of conditions - accretion disks in black hole binaries 
are truncated quite sharply at the ISCO radius. In this section, we first discuss the 
physics of these disks, and we close by describing the mechanics of continuum fitting. 



3.1 Accretion Disk Theory 



The basic physics of black hole acc retion is straightforward ( Frank et al.ll20o3 : lKato et al 



l2008l : lAb ramowicz fc Fragildlioii '). Gas with angular momentum flows in from the out- 
side and settles into a circular orbit stabilized by centrifugal force. The gas steadily loses 
angular momentum as a result of magnetic stresses from the magnetorotational insta- 
bility (Balbus fc Hawlevi,199Sl. who s e effec t is often approximated via the a-viscosity 
prescription of [shakura fc SunvaevI 1 1973h . As the gas loses angular momentum, it 



moves inward, occupying at each instant a circular orbit appropriate to its instanta- 
neous angular momentum. The inward drift continues until the gas reaches the radius 
of the ISCO, -Risco- Inside -RiscOi no stable circular orbits are available and the gas 
falls dynamically into the black hole. 

As described above, the ISCO represents a major transition point in disk physics, 
where gas switches from slow viscous accretion on the outside to inviscid free-fall on 
the inside. The ISCO is thus effectively the inner edge of the disk. Correspondingly, 
information on the linear dimensions of the radius -Risco is imprinted on the emitted 
radiation. Since -Rts oo varies monotonically with the black hole spin parameter a* 



I Bardeen et al.|[l972l ). as illustrated in Figure it is thus possible to measure a* by 



modeling the disk emission. 

The model of choice for this purpose is that described bv lNovikov fc Thornd (|l973l ') , 



hereafter referred to as the N T model, which is the rela tivistic generalization of the 
thin accretion disk model of I Shakura fc SunvaevI (|l973h . Using nothing more than 
the Kerr metric, basic conservation laws of mass, momentum, angular momentum 
and energy, and assumptions of axisymmetry a nd steady state, the NT model (see 
also IPaee fc Thom3l 19741 : iRiffert fc Heroldll 19951 ') derives an analytical formula for the 
differential luminosity dL{R) / dR emitted by the disk as a function of radius R. 

The solid lines in Figure (Sja show for three values of a* the differential disk lumi- 
nosity predicted by the NT model. The disk flux vanishes at -Risco because the model 
has, by assumption, no viscous stress inside this radius (see Section r4.2l for further dis- 
cussion) . More importantly, the peak emission occurs at a radius that tracks the ISCO 
(it is a factor of ^ 2 — 3 larger than -Risco)- This means that the radiation is emitted 
from a progressively smaller effective area, roughly oc -Rjgcc the black hole spin 
increases. Therefore, for a given total disk luminosity, the temperature of the emitted 
radiation increases with increasing a*. This is the key physical effect that underlies the 
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a. R/[GM/c^) 

Fig. 3 (a) Radius of the ISCO Risco s-nd of the horizon ijjj in units of GMjc? plotted as a 
function of the black hole spin parameter a*. Negative values of a* correspond to retrograde 
orbits. Note that -Risco decreases monotonically from 9GM/c^ for a retrograde orbit around 
a maximally spinning black hole, to &GM/c^ for a non-spinning black hole, to GMjc? for a 
prograde orbit around a maximally spinning black hole, (b) Profiles of d{L/ M)/dln R, the 
differential disk luminosity per logarithmic radius interval normalized by the mass accretion 
rate, versus radius R/ (GM/c^) for three values o f a, . Solid lines are the predictions of the NT 
model. The dashed curves from lZhu et al ] 1I2OI2II . which show minor departures from the NT 
model, are discussed in Section 15.21 



continuum-fitting method. By measuring the characteristic temperature and luminos- 
ity of the disk emission, and applying the NT model, one is able to estimate both a* 
and the mass accretion rate M. 

As should be clear from the above, the accuracy of the continuum-fitting method 
ultimately depends on the reliability of the NT model; this issue is discussed further in 
Sections l4.2l andl 5.2l It also depends on our ability to calculate the spectrum of the radi- 
ation, which would be trivial if the disk radiated as a perfect blackbody. Unfortunately, 
because electron scattering plays a prominent role at the X-ray temperatures found in 
black hole binaries, the emitted spectrum is substantially harder than a blackbody 
spectrum of the same flux. Hence it is necessary to employ detailed disk atmosphere 
m odels. Most of the work t o date is based on the atmosphere model bhspec developed 
bv lPavis fc Hubenvl ( 20061 ) . which is discussed in Sections 3.2 and 5.3. 



3.2 Continuum Fitting in Practice 

In broad outline, one fits the X-ray continuum spectrum to the Novikov-Thorne model 
of a thin accretion disk with other spectral components as needed, principally a Comp- 
ton component. As stressed in Sections 1 and 5, in order to obtain useful constraints 
on a*, one must inform the fitting process by inputting accurate value of the external 
parameters D, i and M. The spectral fit returns two output parameters: the spin a* 
and the mass accretion rate M. An important derived quantity is the Eddington-scaled 
luminosity of the disk component L(a*, M)/Z/Edd- 
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In practice, one usually fits the thermal component usi ng; kerrbb2 ( McClintock et al.l 
bood 'n. which is a hybrid code implemented in XSPE C JArnaudl 1996 )that combines 
the capabilitie s of two relativistic disk models, bhspec ( Davis et alj 20051 ') and kerrbb 
( Li et alJliooi ). This latter model, kerrbb, which is a straightforward implementation 
of the analytic Novikov-Thorne model, has three principal fit parameters: a*, M, and 
the spectral hardening factor /, which relates the observed color temperature to the 
effective temperature, / = T/T^g. 

In fitting the disk component with kerrbb, it is quite generally the case that one 
can only determine two parameters, a shape parameter (e.g. a* or T) and a normal- 
ization constant (e.g., M). That is, in practice one cannot additionally obtain a useful 
constraint on /. However, this limitation of kerrbb is handily overcome by pairing 
it with BHSPEC, which is based on non-LTE disk atmosphere models within an a- 
viscosity prescription, bhspec has just two principal fit parameters (spin and mass 
accretion rate) and can be used to fit directly for a*. However, it does not include the 
effects of self-illumination of the disk ("returning radiation"), which is a feature that 
is included in kerrbb. 

The pairing of kerrbb and bhspec is achieved using kerrbb2, which is a modified 
version of kerrbb that contains a pair of look-up tables for / corresponding to two 
values of the viscosity parameter: a — 0.01 and 0.1. The entries in the tables are 
computed using bhspec. The two tables give / versus -£//-£/Edd for ^ wide range of the 
spin parameter (|a*| < 0.9999). The computations of / versus L/I/Edd ^re done using 
the appropriate, corresponding response matrices and energy ranges used in fitting the 
spectra with kerrbb. Thus, kerrbb and the subroutine/table computed using bhspec 
(which together constitute kerrbb2) allow one to fit directly for a* and -£//-£/Edd) while 
retaining the returning-radiation feature of kerrbb. 

Depending on the quality of a particular spectrum, it may be necessary to include 
minor spectral components (e.g., line or edge features), but these cosmetic features 
do not significantly affect the spin results. The three key model components, which 
are fitted in conjunction with the thermal com ponent, are a low-energy cutoff, a "re- 
flected" component (e.g., IRoss fc F abian"2007[l , and a Compton component. The cutoff 
is straightforward to model (e.g., Iwilms ct aU bOOd ). and the reflected component is 
relatively weak i n disk-dominated spectra, even in the most extreme circumstances 
( Gou et aP I2OI1 I). It is the modeling of the Compton component that has been of 
central concern in applying the continuum-fitting method, and we discuss this issue 
now. 

All spectra of black hole binaries, even the most disk dominated, show a high- 
energy tail component of emission, which i s widely attributed to Compto n upscat- 
tering of soft photons by coronal electrons (Rcmillard fc McClintockl bood '). In early 
continuum- fitting work (Sh afee et a l. 2006; McClin tock et al.. .20061 ■ this component 
was modeled unsatisfactorily by adding a power-law component to the spectrum. All 
subseque nt work has used a much-improved empirical model of Comptonization called 
SIMPL (S teiner et al.|[2009bh . This model self-consistently generates the Compton com- 
ponent from the thermal seed spectrum of photons. It allows reliable measurements of 
spin to be obtained even as the fraction of seed photon s /gc th at ar e scattered into the 
power-law component approaches 25% ( Steiner et al.l 2009a, 20 09bh . The use of SIMPL 
in place of the standard power law has doubled the body of useful data for several 



^ Fo r alternatives, see lGierlinski et al.l (l200lh . lKoIehmainen fc Pond (I2OIOI) and lStraub et aO 
j2011l'l . 
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Fig. 4 Model fit to a disk- 
dominated spectrum of LMC 
X-3 obtained using detectors 
aboard the BeppoSAX satel- 
Hte for D = 52 kpc, i = 67° 
and M = 10 Mq llDavis et al.l 
2006). A green solid curve, 
which is difficult to discern 
because it hugs the data, is 
the total model. Also shown is 
the thermal component (red 
long-dashed curve) and the 
Compton component (violet 
short-dashed curve). The re- 
flected component is negligi- 
ble and was not included. The 
orange solid curve shows the 
total model with the effects 
of interstellar absorption re- 
moved. Note that the peak 
Compton flux is only 1% of 
the peak thermal flux. 



sources (e.g., see lSteiner et al.1l201ll . [2012al ). and it has enabled the measurement of 
the spins of bl ack holes whose spectra are pers istently quite str ongly Comptonized such 
as LMC X-1 (|Gou et ahlboogl ) and Cyg X-1 (|Gou et ahlboTlh . 

Successful application of the continuum-fitting method requires the selection of 
spectra that are disk dominated. For transient sources, such spectra are typically ob- 
served during the latter part of an outburst cycle (Figure [JJj) . Figure |4] shows a spec- 
trum with a peak flux in the Compton component that is only 1% of the peak thermal 
flux. For spectra that are this disk-dominated, how one chooses to model the adul- 
terating Compton component is obviously quite unimportant. Meanwhile, there is an 
abundance of spectra of comparable q uality available for se veral sources, i.e., sources 
with /sc ~ 1% (e.g., see Figure 1 in ISteiner et aT]|2009al . and Table 1 in Steiner et 
al. 2011). 

While it is e ssential to select sp ectra that have a substantial thermal component 
(i.e., fsc ^ 25%: 'steiner et al.ll2009ah . it is equally important to select data of moderate 
luminosity, specifically spectra with Eddington-scaled disk luminosities I//I/Edd < 0.3. 
Otherwise, the d isk scale-height grows a nd the thin disk model is invalidated {Sec- 
tions and [5^ iMcClintock et ahlbood ). Fortunately, there is usually an abundance 
of such data because a typical transient source remains for months in a suitable disk- 
dominated state of moderate luminosity during the latter part of its outburst cycle 
(see Figure [2}3). A very wide range of detectors are capable of providing suitable data 
(see example in Section [4. ip . The principal requirements are that the data can be cor- 
rected for dead time, and that the detector have a dozen or more energy channels, an 
appropriate bandwidth, and be well calibrated (Section [571} • 



4 Truncation of the Disk at the ISCO 



We review the large body of observational evidence that there exists a constant inner- 
disk radius in disk-dominated states of black hole binaries. We follow with theoretical 
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evidence, based on GRMHD simulations, that this fixed radius can be identified with 
the radius of the ISCO. 



4.1 Observational Evidence 

It has been clear for decades that fitting the X-ray continuum might prove to be 
a promising approach to measuring black hole spin. The earli est indications came 
with the advent in the mid-1980s of a nonrelativistic disk model ( Mitsuda et al.lll984l : 



Makishima et allllQSfih . now referred to as diskbb, which returns t he color tempera- 
ture T in at the inner-disk radius Ri^. In an important review paper. IXanaka fc LewirJ 
( 19951 ) show the remarkable stability of R^^ for three transients as the thermal fiux of 



these sources steadily decays on a timescale of months by factors of 10-100 (see their 
Figure 3.14). Tanaka & Lewin remark that the constancy of Ri^ suggests that this fit 
parameter is related to the radius of the ISCO. Subsequently, similar evidence for a con- 
stant inner-disk radius in disk-dominated states of black hole binaries has been demon- 
strated for many sources by showing that the bolometric luminos ity of the thermal com- 



ponent is approximately proportional to Tjt ( Kubot^_et_al. |2001; Ku bota fc Makishimal 
I2OO4I : iGierlihski fc Donell2004l : lAbe et alEoOSl : iMcClintock et al... 20091 )! 



A recent study of the persistent source LMC X-3 presents the m ost compelling 



evidence to date for a constant inner-disk radius ( Steiner et aLllioiol l. This result is 
based on an analysis of a large sample of X-ray spectra collected during eight X-ray 
missions that span 26 years. As illustrated in Figure [5] for a selected sample of 391 
RXTE spectra, the radius of the accretion disk was found to be constant over time and 
unafi^ected by the gross variability of the source to within ~ 2 percent. Even considering 
an ensemble of eight X-ray missions, the radius was observed to be stable to within 
« 5 percent. These results provide compelling evidence for the existence of a fixed 
inner-disk radius and establish a firm empirical foundation for the measurement of 
black hole spin. The only reasonable inference is that this radius is closely associated 
with the radius of the ISCO, as we show to be the case in the following section. 



4.2 Theoretical Evidence 

The NT model makes one key assumption: It assumes that the viscous torque vanishes 
inside the ISCO. While em inently reasonable (e.g., .Pa.czvrlski 2000 ; Afshordi fc Paczvhsk| 
l200i IShafee et al.ll2008bl ) , this "zero-torque" assumption does not follow directly from 
basic conservation laws but is applied as an extra ad hoc boundary condition. Further- 
more, the luminosity profiles shown for the NT model in Figure [3]3 depend critically on 
this boundary condition because this condition causes the lumino sity profiles t o vanish 
at the ISCO, which in turn fixes the radius of peak disk emission. iKrohkl l|l999l ) argued 
that magnetic stresses can operate freely across the ISCO an d will cause stron g torques 
at the ISCO, as well as in the inner plunging region. Also, ICammid |l99i) came up 



with a simple analytical MHD model of the plunging region with demonstrably non- 
zero torques. What do real disks do? 

To answer this question, geometrically thin accretion disks around black holes have 
been simulated b y a number of authors (Shafco ot al. 2008a, ; Nobl e & Krolik 20091; 
IPenna et allboioh using state-of-the-art GRMHD codes. The main advantage of simu- 
lations is that they do not require ad hoc assumptions. One simply introduces magne- 
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Fig. 5 {top) Accretion disk luminosity in Eddington-scaled u nits (for M = 10 Mp ,) versus time 
for all the 766 spectra considered in a study of LMC X-3 bv lSteiner et al.1 1120101) . (Downward 
arrows indicate data that are off scale.) Selected data in the unshaded region satisfy the 
thin-disk selection criterion L/LEdd < 0-3 and avoid confusion with stron gly-Comptonized 
hard-state data with fsc > 25% (Section 3.2; lRemillard fc McCnntockll2006l ) . (bottom) Fitted 
values of the inner-disk radius are shown for thin-disk data in the top panel that meet the 
selection criteria of the study (a total of 411 spectra). Despite large variations in luminosity, 
Tin remains constant to within a few percent over time. The median value for just the 391 
selected RXTE spectra is shown as a red dashed line. 



tized gas in a Kerr space-time and lets the system evolve to a quasi-steady state. Since 
the ISCO lies inside the simulation box, well away from computational boundaries, it 
is not treated differently from other regions of the system. In other words, no boundary 
condition is applied by hand at the ISCO, a great improvement over analytical models. 
On the other hand, for technical reasons, simulations to date have not treated radiation 
transfer self-consistently but instead have assumed local cooling. This is not considered 
serious for the purposes of testing the zero-torque condition. 

The dash ed lines in Figure [3b show results from simulations o f very thin disks 
{H/R ~ 0.05:IPenna ct al."201(]l: lKulkarni et~al]l201ll : [zhu et al.ll2012l ). The simulation- 
derived disk luminosity profiles show modest deviations from the NT model predictions; 
in particular, the disk flux does not vanish inside the ISCO. On the other hand, the devi- 
ations are minor, and we discuss and quantify these effects in Section [5]2] Importantly, 
the radius corresponding to the luminosity peak, which is the most relevant quantity 
for the continuum-fitting method, agrees quite well. As discussed in Section 15.21 the 
good agreement between model predictions and simulation results translates into mod- 
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est uncertainties in spin estimates. GRMHD simulation results such as these shown in 
Figure |3}3 are viewed as a strong validation of the NT model. 

Interestingly, t he magnetic s t ress i n the simulations does not vanish in the plung- 
ing region. Indeed, iPenna et al ] (|20ld ) found that the stress there agrees remarkably 
well with Gammie's (1999) model. However, there is little energy dissipation associ- 
ated with this stress (Gammie's analytical model has zero dissipation), so it has little 
bearing on the continuum-fitting method. Another interesting result is that deviations 
from the NT model seem to increase as the luminosity and, concomitantl y, the disk 
thickness increases jKulkarni et all 1201 ll ). as anticipated in previous work (|Paczvriskil 
l200(]| : IShafee et al.ll2008bh . 



5 Uncertainties in Spin Estimates 

The bottom line of this section is that the error in a* is dominated by the observational 
errors in the external input parameters D, i and M. By comparison, the errors due to 
reliance on the NT model, as well as on the disk atmosphere model that delivers the 
spectral hardening, are less important. Meanwhile, the one significant question that 
hangs over most of the spin results is the assumption that the black hole's spin vector 
is aligned with the orbital angular momentum vector. We discuss these points in turn. 



5.1 Observational Errors 



In early work, the error in a* attributable to the uncertainties in D, i and M w as only 
crudely estimated (jShafee et allbood : loTvis et al.ll2006l: iMcClintock et al.|[20oi ). How- 



ever, in all subsequent work, starting with 



Liu et al 



20081 ) , the error in the spin due to 



the combined uncertainties in these three parameters has been computed in detail via 
Monte Carlo simulations. That the error budget for a* is dominated by the uncertain- 
ties in D, i and Al has been thoroughly demonstrated in recent work, which provides 
error estimates for a very wide range of statistical and systematic errors associated 
with (1) the details of the spectral models employed, (2) flux calibration uncertainties, 
(3) the effects of a warm absorber, etc. Rather than discuss these technic al details here, 
w e refer the reader to Section 5 and Appendix A in lsteiner et al ] (|201ll ) and Section 5 
in lGou et all (|201lh . 



5.2 Errors from the Novikov-Thorne Model 

As described in Section [3.11 the NT model is robust and makes very few untested as- 
sumptions. It is true that some properties of the disk, e.g., the density and temperature 
of gas at the disk mid-plane, depend on the magnitude of the viscosity parameter a, 
but the all- important luminosity profiles shown in Figure [HJa do not. These profiles are 
a direct consequence of energy conservation - gas drifts inward, it converts gravita- 
tional potential energy into orbital kinetic and gas thermal energy, and the latter is 
radiated. This physics is independent of the value of a, or even the validity of the a 
prescription. 

The NT model assumes that dissipated energy is radiated by the gas locally at 
the same radius. This is a very safe assumption. The cooling time of the gas in a thin 
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accretion disk is approximately (H/R)'^ times the viscous radial advection time. For 
the disks of interest to the continuum-fitting method {H/R < O.f), this means that 
cooling is about 100 times faster than energy advection and hence very local. 

Another approximation in the NT model is the neglect of disk self-irradiation. 
This is acceptable near the peak of the luminosity profile, where local energy dissi- 
pation greatly exceeds irradiation. However, it is less safe at larger radii. The models 
KERRBB and KERRBB2 (Section I3.2|l include self-irradiation consistently. In practice, 
self-irradiation seems to have a minor effect on spin estimates. 

As already discussed, the NT model assumes zero torque at the ISCO. Although this 
approximation turns out to be less severe than one might have anticipated (Figure Od), 
we still expect it to have some effect on the continuum-fit ting method. Several authors 
(jKulkarni et al.ll201ll : iNoble et ai]l201ll : IZhu et al.ll2012h have investigated this issue 
quantitatively. The general consensus is that the zero-torque approximation introduces 
uncertainties in spin estimates of around Zia* 0.1 for low spin values a* < 0.5 and 
much smaller errors as a* 1. These error estimates are for geometrically thin disks: 
H/R < 0.1, or I//I/Edd ^ 0.3. The errors are more severe for thicker disks. 

As a bottom line, these estimates of the NT model errors are significantly less 
than the observational errors presented in Section |6l Furthermore, these NT model 
errors are overestimates because most the spin measurements to date are based on 
spectral data with L/Lgj^j < 0.1, and, correspondingly, data for disks so thin that the 
MRI turbulence cannot now be resolved in practical GRMHD simulations. Looking to 
the future, it might be possible to do better by replacin g the NT model with a more 
accurate simulation-based model (e.g.. |Penna et ani2012h . but this step is not presently 
warranted. 



5.3 Errors from the Disk Atmosphere Model 

An essential cornerstone of the contiimum-fit ting method is a reliab le model of the 



disk's atmosphere. Such a model is bhspec jPavis fc Hubenvlbood ). which can be 
used either alone or including the effect of self-irradia tion via ker rbb2. bhspec is 
based on the non-LTE radiative transfer code tlusty (|Hubenv fc L anz 1995), which 
was originally developed for stellar atmospheres. It is quite sophisticated and includes 
a wide range of physical effects. 

At a given location on an accretion disk, bhspec computes the emitted spectrum 
using three supplied parameters: the effective temperature T^g defined such that ra- 
diative fiux F — crT^ff, local vertical gravity parameter Q, and disk column density 
S. As discussed in previous sections, a robust estimate of T^g can be obtained from 
the NT model, while the parameter Q is calculated directly from the Kerr metric. The 
main uncertainty is in the value of S. 

In standard disk theory, S varies inversely as the viscosity parameter a and is thus 
quite uncertain. Fortunately, in the case of optically thick disks (which all thermal state 
disks are) S has only a weak effect on the emerging spectrum. This is analogous to the 
case of a star where the spectrum depends on the effective temperature and surface 
gravity, but not at all on the optical depth to the stellar core, which is effectively infinite. 
The optical depth through a disk is not quite infinite, hence there is some spectral 
dependence on S. However, this dependenc e is weak for models with L/LEdd < 0-3 
( Davis fc Hubenvlbood : [Pone fc Davidboosl ). 



14 



For the same reason, details of exactly how viscous heating is distributed vertically 
within the disk are unimportant. So long as energy dissipation occurs in the disk 
interior at o ptical depths greater tha n a few, the emerging spectrum depends only on 
T^ff and Q I Davis et al-lbOOsl . l2009l l . This is not true if there is substantial energy 
dissipation close to or above the photosphere. Disks in the thermal state probably do 
not have such dissipation since their spectra show very little hard "coronal" emission 
( Remillard fc McClintoc5 '2006''). Whatever little coronal emission is present is fitted 
for via a model for the Compton power law such as SIMPL (Section 13. 2p . 

The standard bhspec model assumes hydrostatic equilibrium a nd does not include 
the force from magnetic fields. However, numerical simulations (e.g.. lHirose et aPboOSh 
indicate that the photospheric surface regions show modest deviations from hydrostatic 
equilibrium and are primarily supported by magnetic forces. Including these effects 
i n BHSPEC genera lly leads to a modest (< 10%) increase in the spectral hardness 
( Davis et al.l 20091 '). Meanwhile, the effects of irradiation (both self- and from a corona), 



which have not yet been rigorously explored, may lead to a slight hardening of the 
spectrum. In summary, while there are uncertainties associated with the disk spectral 
model used in the continuum-fitting method, it appears unlikely that the resulting 
errors in J^in are more than 10%, which for low values of spin implies Z\a* ~ 0.1, 
decreasing as a* ^ 1 (Figure [3^). 



5.4 Assumption of Spin-Orbit Alignment 

In determining the spins of eight of the ten black holes (see Section |6}, it is assumed 
that the plane of the inner X-ray-emitting portion of the disk is aligned with the binary 
orbital plane, who se inclination angle i is determined from optical observations (e.g., 
lOrosz et al.ll2011al ). However, if a black hole's spin is misaligned with the orbital vector, 
this will warp a thin disk because the Bar deen-Petterson effect will force the inne r disk 
to align with the black hole spin vector ( Fragile et al.l[2007l : iMartin et aPboOTl 'Fl An 



error in estimating the inclination of the inner disk of ~ 10° or more, resulting from an 
erroneous use of i as a proxy for the inclination of the inner disk, would substantially 
corrupt most continuum-fitting measurements of spiiJ^J. 

There is evidence for gross spin-orbit misalignm ent for one t ransient system (SAX 
J1819. 3-2525); however, this evidence is weak (iNaravan fc McCl intock 2005). F or the 
transients generally, more recent evidence, which is summarized in Section 1 of lsteiner fc McClintockl 
( 20121 ). argues in favor of alignment. Briefiy, the timescale for accretion to torque the 



black hole into alignment is estimated to be ~ 10^ — 10^ years, which is short com- 
pared to the typical lifetime of a transient system (Section 2.2). In the case of the 
persistent supergiant systems, there is some evidence that their more massive black 
holes are formed by direct, kickless collapse (|Mirabel fc RodriguesI [20031 : | Reid et al.l 
l201lh . Finally, a population synthesis study based on a maximally conservative (i.e., 
minimum-torque) assumpti on indicates that th e spin axes of most black hole primaries 
will be tilted less than 10° (jFragos et al.ll2010l) . 



In determining the spins of the remaining two black holes (see Section [S}, the 
inclination of the inner disk is taken to be the inclination of the radio or X-ray jet axis, 



* While thin disks are subject to warping, thick disks are not (iDexter &: Fragi"i3l201ll V 
^ Unfortunately the continuum-fitting method cannot fit for the inchnatio n of the inner d isk 
because there is a degeneracy between the inclination and spin parameters jLi et al.ll2009f) . 
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which is presumed to be aligned with the black hole's spin axis. The jet inclination angle 
for these microquasars, GRS 1915+105 and H1743-322, was de termined by modeling 
prop e r-motion data derived from radio and X -ray observations I Mirabel fc Rodriguea 
ll994l : lFender et al.|[l999l : ISteiner et al.ll2012al '). Fortunately, radio/X -ray jet data have 
also yielded a strong constraint on 9 for a third microquasar, XTE J1550-564, thereby 
providing a rare opportunity to check directly the assumption of spin orbit alignm ent 
because i ts orbital inclinat ion angle i has also been measured ( Orosz et al.ll20l"lbh . In 
this case, ISteiner fc McClintock (.2012 ) find no evidence for misalignment and place an 
upper limit on the difference between the spin and orbital inclinations of \0 — i\ < 12 deg 
(90% confidence). 



6 Results and Discussion 

Table 1 lists the masses and spins of ten stellar black holes. By virtue of the no-hair 
theorem, this table provides complete descriptions of each of these ten black holes. 
The spins span the full range of prograde values, and the masses range from 6 to 
16 M p,. In addition to the continuum- fitting spin data in Table 1, ICierlinski et al.l 
provide prelim i nary e stimates for the spins of LMC X-1 and GRO J1655-50, 



(12001 



iKolehmain en fc Done (20W) report a hard upper limit of a* < 0.9 on the spin of 
GX 339-4, and Nowak ct al. (2012|) argue that the spin of 4U 1957-1-11 is extreme. 
Concerning 4U 1957-1-11, it is unclear if the compact object i s a black hole, and the 
key parameters D and M are essentially unconstrained. Finally. [Middleton et al. I diooi ) 



find an apparently moderate value of spin for GRS 1915-1-105, which is at odds with 
the extreme value in Table 1; the reason why Middleton et al. ob tained a depressed 
value of spin is explained in Section 5.3 in lMcClintock et al.l 1 2006h . 



6.1 The Persistent Systems vs. the Transients 

Ignoring for now GRS 1915+105, there is a striking dichotomy between the black 
holes in persistent systems and those in transients, both in their masses and their 
spins (Table 1). Considering spin first, the three persistent black holes all have high 
or extreme spins. In contrast, the spins of the six transient black holes (i.e., excluding 
GRS 1915+105) range downward from ~ 0.8, while the spins of four of them are consis- 
tent with zero. The dichotomy is even sharper if one considers the spins predicted for six 
additional transient black holes, which are based on a fitted correlation between radio 
powe r and spin (Section l7.2p : their spins are all predicted to be a* < 0.8 (jSteiner et al.l 
I2OI3I V Concerning the extreme spin of the transient GRS 1915+105, we simply note 
that this prototype of the microquasars is a remarkable transient for several reasons, 
including its relatively high mass, unique X-ray variability, constellation of four high- 
frequency X-ray oscillations, extreme 33.9-day orbi tal period, and maintenance o f a 
high X-ray luminosity since its appearance in 1992 (|Remillard fc McClintoc5l2006l ). 



Not only are the persistent black holes rapidly spinning, they are also massive - 
11 — 16 Mq ~ compared to the transient black holes. The masses of the transients are 
significantly lower and, remarkab ly, their mass distribution is narrow: 7.8 ± 1.2 Mq 
( 6zel et alllioiol : iFarr et al.ll201lh . 



16 



Table 1 The masses and spins, measured via continuum- fitting, of ten stellar black holes". 



System 




M/Mq 


References 


Persistent 


Cyg X-1 


> 0.95 


14.8 ± 1.0 


Gou et al. 2011: Orosz et al. 2011a 


LMC X-1 




^■^^-0.07 


10.9 ± 1.4 


Gou et al. 2009; Orosz et al. 2009 


M33 X-7 


0.84 ±0.05 


15.65 ± 1.45 


Liu et al. 2008: Orosz et al. 2007 


Transient 


GRS 1915+105 


> 0.95'' 


10.1 ± 0.6 


McClintock et al. 2006; Steeghs et al. 2013 


4U 1543-47 


0.80 ±0.10'' 


9.4 ± 1.0 


Shafee et al. 2006; Orosz 2003 


GRO J1655-40 


0.70 ±0.10'' 


6.3 ± 0.5 


Shafee et al. 2006; Greene et al. 2001 


XTE J1550-564 


34+0-20 


9.1 ± 0.6 


Steiner et al. 2011: Orosz et al. 2011b 


H1743-322 


0.2 ± 0.3 


~ 8^= 


Steiner et al. 2012a 


LMC X-3 


< O.S'* 


7.6 ± 1.6 


Davis et al. 2006: Orosz 2003 


A0620-00 


0.12 ±0.19 


6.6 ±0.25 


Gou et al. 2010: Cantrell et al. 2010 


Notes: 

" Errors are quoted at the 68% level of confidence, 
estimated to be at the 99% level of confidence. 


except for the three spin limits, which are 



'' Uncertainties greater than those in papers cited beca use early error es timates were crude. 
" Mass estimated using an empirical mass distribution llOzel et al.ll2010h . 
Preliminary result pending improved measurements of M and i. 



6.2 Prograde Spins that obey the Kerr Bound 

The lack of negative spins in Table 1 may be the result, in a close binary system, of the 
expected alignment of the spin of the black hole progenitor with the orbital angular 
momentum, and it may also indicate that black hole kicks are not strong enough to 
flip the black hole into a retrograde configuration. While interesting that there are no 
negative spins, it is equally interesting that the spins of all ten black holes obey the 
Kerr bound a* < 1. In particular, if the distances to either Cyg X-1 or GRS 1915+105 
were ~ 30% less than the best current estimates, then it would be impossible to fit the 
data with the kerrbb2 model, which only accommodates spin values a* < 1. Because 
the observed values of each of the three external fit parameters {D, i and M) place 
hard constraints when fitting the data, a failure to fit a spectrum that requires a* > 1 
h as the potential to falsify the spin model. For a discussion of this point, see Section 6.4 
in lMcChntock et al.l (|2006h . 



6.3 The High Natal Spins of the Persistent Black Holes 

It is reasonable to conclude that the black holes in the persistent systems were born 
with high spins because their host systems are too young for these black holes to have 
be en spun up by ac cretion torques. Consider, for example, the persistent system Cyg X- 
1 (|Gou et al.ll20Til : For its black hole to achieve its present spin of a, > 0.95 via disk 
accretion, an i nitially nonspinning black hole wo uld have had to accrete > 7.3 Mq 
from its donor (|Bardeenlll970l : lKing fc Kolblll999l ) to become the 14.8 Mq black hole 
we observe today. However, even at the maximum (Eddington-limited) accretion rate 
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this would re quire > 31 million y ears, while the age of the system is between 4.8 and 7.6 
million years (jWong et al.ll2012l l. Likewise for M33 X-7 and LMC X-1, the correspond- 
ing minimum spin-up timescales are > 17 and > 25 million years, respectively, whil e 
the respective ages of the systems are < 3 and < 5 million years ( Gou et al] 201 il l. 
It therefore appears that the spins of these systems must be chiefly natal, although 
possibly such high spin s could be achieved during a short-lived evolutionary phase of 
hypercritical accretion ({Moreno Mendez et a l. 2008). A similar argument for natal spin 
can be made for black holes in transient systems, especially those with relatively high 
spins and massive secondaries, but it is less clear-cut. 



6.4 Applications 

The data in Table 1 have a number of applications to physics and astrophysics, both 
immediate and potential. In physics, a high goal is to use such data as a springboard 
to test the no-hair theorem (see Section [8|, and the foundation for any such test 
is high-quality measurements of mass and spin for a good sample of black holes. In 
astrophysics, knowledge of the spins of stellar b lack holes is crucial for example in con- 
straining models of gamma-r ay burst sources llWooslevlEiii : iMacFadven fc WoQslevI 



1999'; 'Wooslev fc Hegeilbood ): supernovae and black hole formation (Lee et al."2002'; 
Wong ct al. 20 13 ): exotic black hole states and state transitions ( Rcmillard & McClintoc^ 
20061 ): a nd in informing LIGQ /V IRGO modelers who are computing gravitational- wave 



signals ( Campanelli et aT]|2006l l. A central question, which we turn to in the next sec- 



tion, is the role of spin in powering jets. 



7 Jet Power and Black Hole Spin 

Since the spin parameter a* is one of only two numbers that completely characterize 
a black hole (mass M being the other), it stands to reason that it should influence at 
least some observational properties of the hole. The most widely discussed connection 
is to relativistic jets. 

The story goes back to lPenros3 (| 19691 ) who showed that a spinning black hole has 
free energy that can in principle be tapped by specially prepared infalling particles. 
Although Penrose's specific proposal is not considered promising, the idea of extracting 
energy from spin n ing b lack hole s has stuck and has beco me popular in astrophysics . 
iRuffini fc WilsonI (| 19751 . see also lOamour et al.lll978l ) and Isiandford fc Znaied (|l977l ) 



suggested a specific mechanism whereby a force-free poloidal magnetic field around a 
spinning black hole is twisted by frame dragging, thereby producing outgoing Poynting 
fiux along twin jets. We refer to this as the generalized Penrose process. 

GRMHD simulations of accreting black h oles have found MHD jets forming spon- 



taneously from generic initial conditions (e.g ., Koide et al.ll2002l: McKi^nne y fc Gammid 
20041 : iMcKinnevllioOsI : [ Beckwith et al.ll2008l : iMcKinnev fc Blandfordll2009l ). Moreover, 



in one particular simulation involving a rapidly spinning black hole and a strong 
poloidal field, Tchckhovskoy ct al. (2011) showed that the power carried by the jet 
exceeded the total rest mass energy of accreted gas, meaning that the jet extracted 
energy from the spinning black hole. 
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On the observational front, until recently there was no empirical evidence for a 
connection between black hole spin and relativistic jets. We discuss here the first such 
evidence. 



7.1 Two Kinds of Jets in Black Hole Binaries 



Fender et al.l ( 20041 ') identified a number of systematic properties in the radio emission 



of black hole binary jets. They showed that there are two kinds of jets, which we refer 
to as "steady jets" and "ballistic jets," each associated with a specific spectral state of 
the X-ray source. Although we discuss both kinds of jets for completeness, our focus 
here is on the ballistic jet. 

The steady jet is observed as a continuous outflow of plasma in the hard spectral 
state (for a discussion of spectral states in black hole binaries, see iRemillard &: McClintockl 

This jet is small-scale, being observable only out to a few tens of AU, and it ap- 
pears not to be very relativistic. It is present at the very start of a transient's outburst 
cycle. Referring to the X-ray light curves for XTE J1859-I-226 in Figure [2^, the jet is 
present during the first few days when the hard flux (panel a) is most intense and the 
2-12 keV flux (panel b) is increasing rapidly. It then disappears, and it returns only 
near the end of the outburst cycle (beyond the right edge of the plot). The steady jet 
is seen in all transients at low values of M. 

The far more dr amatic ballistic jet is launched when a transient goes into outburst 
(|Fender et al.ll2004 ). This powerful transient jet usually appears near (or soon after) 
the time of outburst maximum, as the source switches from its initial hard state to a 
soft state via the "steep power-law" state. Ballistic jets manifest themselves as blobs 
of radio (and occasionally X-ray) emitting plasma that move ballistically outward at 
relativistic speeds (Lorentz factor F > 2). They are often observed out to distances of 
order a parsec. Because ballistic jets resemble the kpc-sc ale jets seen in quasars, blac k 
hole binaries that produce them are called microquasars (^ Mirabel fc Rodn'guej|l999l '). 

Ballistic j et ejection occurs at a very speciflc stage during the spectral evolution of a 
given system ( Fender et al.l[2004l '). In Figure [J] the strong spike in the radio light curve 



(panel d), which is characteristically delayed relative to the corresponding spike in the 
X-ray luminosity (panel b), is associated with a ballistic jet. As most clearly de mon- 



strated for the prototypical microquasar GRS 1915-fl05 (jFender fc Bellonill2004l l. this 
ejection stage appears to correspond to the inward-moving inner edge of the accretion 
disc reaching the ISCO, which apparently results in some violent event that launches 
a large-scale relativistic blob. 

On general principles, one expects jet power to depend on a black hole's mass M 
and spin a*, and the mass accretion rate M, plus other factors such as the strength 
and topology of the magnetic fleld. If one wishes to investigate the dependence of jet 
power on a*, one needs flrst to eliminate the other variables. 

For steady jets, M spans a wide range, and it is not straightforward to eliminate the 
effects of this variable. It is p ossible to do this, however, by using the disk X-ray lumi- 
nosity as a proxy for M (e.g., Heinz_&_Sunyae^ 20031 : 1 Merloni et al.|[2003l : iFalcke et al.l 
I2OO4I : iFender et al.ll20oi : lFender fc Bellonill2004h . but one must have knowledge of the 



radiative efficiency of the disk, which is general ly both low and v ariable in the hard 
state where the steady jet forms (see lNaravan fc McClintockl l2008l . for a discussion of 
radiatively inefficient accretion in the hard state). The pro cedure is relat i vely u ncer- 
tain and it is difficult to obtain robust results. Nevertheless, iFender et al.l l|2004 ) have 
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performed such a study and have claimed that there is no evidence in the data that 
the power of a steady jet depends on spin. 

Bal Ustic jets on the oth er hand invariably occur near the peak of transient out- 
bursts. ISteiner et al] l|2013h have shown that during major outbursts the peak disk 
luminosities in various transients are near the Eddington limit and are clustered within 
a factor of ~ 2 in luminosity, which means that these systems behave for all purposes 
like "standard candles." This crucially allows one to compare the power of ballistic jets 
observed for different black holes at the same M, namely M ~ Af^dd- In a ddition, all 
black holes in transients have similar masses to better than a factor of two ( Ozel et al.l 
furthermore, it is easy to correct for this small effect (see below). This leaves a* 
(with magnetic field as a wild card) as the sole remaining parameter that could have 
any influence on jet power. 

Let us define the jet efficiency factor 77 of a ballistic jet, 

VjetM = (ijct)/{M)c', (1) 

where (I/jct) is the time-average kinetic luminosity flowing out through the jet and 
{M)(? is the time-average rate at which rest- mass energy flows into the bl ack hole. 
Using (1) radio luminosity as a proxy for Ljet jNaravan &: McClintocQl2012l ) and (2) 
observed va lues for the peak rad io luminosities for five black holes that are all accreting 
at ~i\fEdd jSteiner et ahllioiih . one can infer directly how jet luminosity depends on 
spin, as we discuss in the following section. 



7.2 Correlation Between Spin and Ballistic Jet Power 

A typical ballistic jet blob is initially optically thick and has a low radio power. As 
the blob moves out and expands, the larger surface area causes its radio power to 
increase. This continues until the blob becomes optically thin, after which the flux 
decl ines rapidly. The ov e rall behavior is generally con sistent with an expanding conical 



jet (|van der LaarJ[T966l : iHiellming fc JohnstonI 19881 ). Moreover, as discussed in Sec 



tion 17.41 the peak radio luminosity is expected to scale more or less linearly with the 
jet kinetic energy or kinetic luminosity. Thus, peak radio luminosity is a good proxy 

for jet kinetic luminosity. 

Naravan fc McC lintocic ('2012') considered the peak radio luminosities of ballistic jet 



blobs in four transients, A0620-00, XTE J1550-564, GRO J1655-40, GRS 1915+105, 
and showed that they correlated well with t he correspond i ng bla ck hole spins measured 
via the continuum-fitting methoc0. Later, ISteiner et al.l (|2013l ) included a fifth tran- 
sient, H1743-322, whose spin had been just measured. Figure |6^ shows a plot of the 
black hole spins of these five objects versus a measured quantity called "Jet Power," 
which refers to the radio luminosity vLy = {vSu)D'^ /M (here, not corrected for beam- 
ing), where v = 5 GHz is the radio frequency. Si/ is the flux density in Jy at the peak 
of the ballistic-jet radio light curve, D is the distance in kpc, and M is the black hole 
mass in solar unitfl. That is, the proxy adopted for jet kinetic luminosity is simply 

^ In the case of a fifth transient, 4U1543-47, radio observations did not include the peak of 
the light curve, so one could only deduce a lower limit to the jet power. Note that the radio 
peak can be very narrow in time, e.g., fa 1-day in the case of XTE J1859+226 (Figure[2]l, so 
one requires dense radio monitoring to catch the peak. 

^ The scaling by mass allows for the fact that all the sources are close to the Eddington 
luminosity limit, which is proportional to mass. However, since the masses of the black holes 
differ little (Table 1), the results would be virtually identical if the mass scaling were eliminated. 
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log a. ^isco ' ^ 

Fig. 6 (a) Plot of the quantity Jet Power, which measures the 5 GHz radio luminosity at 
light curve ma ximum, versus black hole spin , measured via the c ontinuum-fitting method for 
five transients l lNaravan &: McClintoclj 120121 : [Steiner et al.|[20T3) . The dashed line has slope 
equal to 2. (b) Plot of Jet Power versus Risco /{GM/c^). Here the radio luminosity has been 
corrected for beaming assuming a bulk Lorentz factor F = 2 (filled circles) or F = 5 (open 
circles). The solid lines cor respond to Jet Power oc O^, where On is the angular frequency of 
the horizon llSteiner et al.ll2013f) . 



the peak radio luminosity at 5 GH j^- Figure shows unmistakable evidence for a 
strong correlation between Jet Power and a*. Note that Jet Power varies by nearly 
three orders of magnitude as the spin parameter varies from « 0.1 — 1. 

The uncertainty in the estimated values of Jet Powe r, which is difficult to assess, i s 
arbitrarily and uniformly assumed to be a factor of two ( Naravan fc McClintockir2012l ) . 
The very unequal horizontal error bars in Figure are a feature of the continuum- 
fitting method of measuring a* . Recall that the method in effect measures -Risco ^nd 
then deduces the value of a* using the mapping shown in Figure [3^. Since the mapping 
is highly non- linear, especially as a* — > 1, comparable errors in -RiscO correspond to 
vastly different uncertainties in a*. In addition, the use of log a* along the horizontal 
axis tends to stretch error bars e xcessively for low spi n values. This point is clarified by 
considering Figure [6j3, based on Isteiner et al.l (|2013t ). Here the horizontal axis tracks 
logi?igco rather than logo*, and the horizontal error bars are therefore more nearly 
equal. The key point is, regardless of how one plots the data, the correlation between 
Jet Power and black hole spin appears to be strong. 



7.3 What Does it Mean? 

Assuming the correlation shown in Figure [6] is real, there are two immediate impli- 
cations: (i) Ballistic jets in black hole binaries are highly sensitive to the spins of 
their underlying black holes, presumably bec ause these jet s derive their power di- 
rectly from the spin ener gy of the hole (a la IPenrosel Il969l : iRufflni fc WilsonI Il975l : 
iBlandford fc Znaieklll977l ). (ii) Spin estimates of stellar black holes obtained via the 
continuum-fitting method are sufficiently reliable to reveal this long-sought connection 
between relativistic jets and black hole spin. 

* None of the results change if one chooses a different reference frequency, e.g., 1.4 GHz or 
15 GHz. 
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With respect to item (i) , we note that the mere existence of a correlation does not 
necessarily imply that the Penrose process is at work. We know that the gravitational 
potential in the inner disk deepens with increasing black hole spin, since the inner 
radius of the disk (-Risco) shrinks with increasing a* (Figure 

Could this disk-related effect be the reason for the increasing jet power? It seems 
unlikely. The radiative efficiency of a Novikov-Thorne thin accretion disk increases only 
modestly with spin; for the five objects shown in Figure |B] the radiative efSciencies are 
0.061, 0.069, 0.072, 0.10 and 0.19, respectively, varying by a factor of ~ 3. It seems 
implausible that a disk-powered jet could vary in radio Jet Power by the orders of 
magnitude seen in Figure [S] 

In contrast, any mechanism that taps directly into the black hole spin energy via 
something like the Penrose process can easily account for the observed large variation 
in Jet Power. For instance, the g eneralized Penrose effect predicts that the jet efficiency 



factor should vary as T^j^t cx: ( Ruffin^_&_^Vilson 1 97 5| : Blandford fc Znaiek|[l977i ) or 



more precisely, as jyjet oc il^ fTche khovskov et al. r i2010l ). where J7h is the angular 
frequency of the black hole horizor[f|. 



^H^^f^^V (2) 



1 + /T 



The dashed line in Figure [6^ corresponds to Jet Power cx a* and the solid lines in 
Figure [61d to cx O^. The observational data agree remarkably well with these scalings, 
strongly suggesting that a Penrose-like process is in operation. 

While the above conclusions are highly satisfying, one should not discount other 
possibilities. First, note that the correlation shown in Figure [6] is based on only five 
objects. Although this is mitigated by the very wide range of Jet Powers, the correlation 
might become weaker with the next spin measurement. Second, the correlation might 
arise if Jet Power and spin are each correlated with some third parameter. For instance, 
it is intriguing that the binary orbital periods of the five transients under consideration 
increase systematically with Jet Power. One could imagine scenarios in which the 
energy of ejected blobs depends on the size of the accretion disk, which depends in turn 
on the orbital period. However, it is less easy to see why the values of a* measured with 
the continuum- fitting method shoul d correlate with the orbital p eriod, while conspiring 
to produce the scaling predicted bv lBlandford fc Znaiekl (|l977l l. 

Returning to the subject of steady jets, it is interesting to co nsider why there is 
apparently no corr elation between jet radio luminosity and spin l|Fender et al.ll2010l : 



iRussell et al.|[20l3 ). One likely answer is that the two kinds of jet are produced via very 
different mechanisms. Perhaps ballistic jets are launched within a few gravitational 
radii, near -RiscOi where the black hole spin could plausibly have a strong effect, 
whereas steady j ets in the hard sta te originate much further out at radii ~ 10 — 
100 GM/c^ ^.i^.- lMarkoff et al.ll2005l ). where the effects of spin are relatively weak. In 
support of this e xplanation, ballistic jets a re definitely relativistic, with Lorentz factors 



of up to several ( Fender et al.l 20041. 20061). whereas there is little ev idence that steady 
jets are relativistic (|Gallo et al.ll2003l : iNaravan fc McClintoc5l2005l '). 



Tile two scalings agree for small values of a, , but differ as a» 
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7.4 A Challenge 

Based on an analysis of a heterogeneous data sample of uneven quality, I Fender et alj 
claim that there is no evidence for a correlation between the power of ballistic 
jets and black hol e spin. The substantial d ifferen ce between the results obtained by 
these authors and iNaravan fc McChntockl (|2012l ) is, in the end, determined by the 
quantity used to represent jet power (which is discussed further in Section I7.4.3[) . 
Fender et. al compute jet power from the peak radio luminosity and rise time of a 
particular synchrotron event; adopt a formula relating jet power to X-ray luminosity, 
namely, logio-f/jct = c + 0.5(logiol/x — 34); and use the normalization constant 
c as their proxy for the jet power. Narayan and McClintock, on the other hand, use 
the model-independent proxy discussed above, namely, the maximum observed radio 

luminosity at 5 GHz. 

Very recently, th e three authors of I Fender et al.l ( 20ld ) have written a second paper 



( Russell et al]|2013l ) repeating their claim that there is no evidence for a correlation 



between the power of ballisti c jets and black hole spin. T herei n, they challenge th e 
methodologies and findings of lNaravan fc McChntocM (|2012l ') and lSteiner et al.l (|2013l ') 



A response to this challenge is being readied (J. Steiner et al., in preparation); mean- 
time, the following is a preliminary sketch of the elements of this response. The foUow- 
ing comments perta in only to ballistic jets and continuum-fitting spin measurements. 
( Russell et al]|2013l . also discuss steady jets and Fe-line spin measurements). 



7.4-1 Significance of the Result 

iRussell et al.l (|2013h contend that the empirical correlation shown in Figure [6^ is only 
marginally statistically significa nt (~ 90% c onfidence). Their analysis is based on a 
Bayesian line a r regr ession model (|Kellvll2007l '). Adopting their procedure, but u sing the 



Steiner et all \20ii ) data set, we confirm the correlation at « 95% confidenc In 
contrast, a traditional analysis (i.e., a linear fit using two-dimensional error bars) pro- 
duces a correlation at ~ 99.9% confidence. 

Of greater importance, iRussell et al.l (|2013h consider only the empirical correla- 
tion shown in Figure [5^, and they wholly disregard the physical model which takes 
into account beaming effects, namely, the model shown in Figure (Gja. The application 
of this physical model is ba sed on four simple assumptions: (i) spin-orbit alignment 
(|Steiner fc McClintockll2012l) : (u) F is the same for all five sources and bracketed be- 



tween 2 and 5 (Fender' 2006^; (iii) jet power is proportional to (|Blandford fc Znaiekl 



[l977. : iTchekhovskov et a l. 2010); and (iv) radio luminosity ca n be u sed as a proxy for 
jet power (Sections 17.21 and 17.4. 3|) . Fitting the ISteiner et all (|2013l ') data set to this 
model, with the model normalization as its sole fit parameter, one obtains good fits 
(shown in Figure [B]d) with /'^ = 0.3 and 0.5 for F = 2 and 5, respectively. This 
relationship is determined over a span of ~ 3 orders of magnitude in jet power, and 
over the full allowed range of prograde spins. It is reasonable to view this result as 
strong evidence that the Blandford-Znajek model successfully describes the behavior 
of a ballistic jet produced by a black hole transient as it approaches its Eddington 
limit. 



Tile lower level of confidence reported by IRussell et al.l l l2013fl is largely attributable to 
their use of radio data for a flare of H1743— 322 tliat is unrelated to the ejection of its ballistic 
jet. 
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Issue of Data Selection 



iRussell et all (|2013l ) furthermore argue that lNaravan fc McClintoc3 (|2012l ') and lSteiner et al. 



( 20131 ) omit data for several systems that should be included in the cor relation plots 



shown in Figure[6l These additional data, which are plotted in Figure Ic in lRussell et alJ 
(|2013l ). destroy the clean correlation shown in Figure |6^. However, it is inappropriate 
to include these data, and we reject them for the following reasons. 

Cygnus X-1 radi ates persistently at a few percent of Eddington and, during its 
periods of jet ejection ( Fender et al .11200^ '). its mass accretion rate is both very low and 



poorly constrained. At the same time, as stressed in Section 17.11 comparing sources 
at a fixed luminosity (namely L ~ L'Edd) is the essential methodological requirement 
that eliminates the otherwise unknown depende ncy of the jet e fBcien cy lyjct on M (see 
equation 1). In order to include Cyg X-1 in the lSteiner et al] (|2013l l sample, it would 
be necessary to know the precise scaling of jet power with M and then to estimate M 
at the time of je t ejection; M i s ~ 1 — 2 orders of magnitude below Eddington and 
highly uncertain. I Russell et al 

] (|2013h ignore this problem. 

For GRS 1124-68 and GS 2000+25, iRussell et aP (|2013l ') adopt continuum- 



fitting spin data that are completely unrelia ble. These data ar e from a pioneering, 
proof-of-concept paper on continuum fitting ( Zhang et al.l IiQQtI ) whose authors note 



that the crucial "system parameters [i.e., D, i and M] are mostly unknown." More 
importantly, the alleged spin estimates for these two sources were simply inferred by 
Zhang et al. from crude estim ates of the inner-disk radius taken from a review 
paper (jTanaka fc Lewinlll995h . and these estimates were computed using the nonrela- 



tivistic disk model diskb b that assumes a grossly incorrect inner-boundary condition 
( Zimmerman et aPboOSl '). while neglecting the effects of spectral hardening. 



GRS J1655-40 had major outbursts in both 1994 and 2005. IRussell et al.1 (|2013h 
plot in their Fig ure Ic a data point for the 200 5 ou tburst which c orresp onds to a very 
faint radio flare. iNaravan fc McClintockl 1 20121 ') and lSteiner et all (2013) do not include 
data for the 2005 outburst in their sample because the radio coverage was too sparse: 
The proximat e observations that bra cket the peak 4.3 Crab flare (see RXTE/ASM plot 
in Figure 1 in lBrocksopp et al]l2006l ) occur 2.3 days before and 4.7 days after the 2-12 
keV maximum (see observation log in the NRAO VLA archive and web link to a plot 
of the radio data in .Rupen et a l. 2005). A bright and brief radio flare, for example like 
the 1-day radio spike observed for XTE J 1859-1-226 (Figure [2Ji), would obviously have 
been missed during this week-long gap in cover age. The radio coverage for t he 2002 



outburst of 4U 1543—47 is likewise inadequate ( Naravan fc McClintockll2012l 'l. 
7.4-3 Synchrotron Bubble Model 

Another point of contention is how to relate the radio luminosity vLu observed at the 
peak of the light cur ve — this is the qua ntity "Jet Power" in Figure 6 — to physical 
conditions in the iet. Isteiner et all ( 20131 ) used synchrotron theory with the following 



standard assumptions: (1) The non-thermal radio-emitting electrons in the jet blob 
have an energy distribution A'^(7) ~ where 7 is the electron Lorentz factor in 
the blob frame; (2) p = 5/2 to be consistent with the synchrotron spectrum; (3) the 
magnetic energy in the blob is in equipartition with the energy in the non-thermal 
electrons; (4) there is one proton for each non-thermal electron and the total energy 
of the blob is dominated by the kinetic energy E of the protons; and (5) at light curv e 
maximum, the jet blob transitions from optically thick to thin ( van der LaanI [T966I ) . 



24 



Steiner et al. showed that E ~ (uLu)^''^, i.e., the blob energy varies approximately 
linearly with vL^. They thus argued that the latter is a good proxy for the former. 

W hy consider blob energy E? In the synchrotron bubble model ( van der LaanI 
Il966l l. the jet ejection is some brief episode that is no t observationa l ly reso lved. Hence 
the total ejected energy is all that one can measure. iRussell et al.l (|2013h (and refer- 
ences therein) focus on the jet kinetic luminosity, Ljet = dE/dt ~ E/tjg^, where tjot 
is the duration of the jet ejection. They further assume that the jet is ejected con- 
tinuously, with a constant lum i nosity , until light curve maximurn^^l. How does Ljot 
depend on vL^? ISteiner et al.l 1 2013h find that the radius of the jet blob scales as 
R ^ {uLu)^'^. Therefore, if the blob expands with some constant speed, say c, then 
Ljet ~ E/{R/c) ~ {vLv)^'"^ , i.e., jet kinetic luminosity varies as vLv to a power some- 
what less than unity. The truth is probably somewhere in between this result and 
that obtained for th e blob energy E, i.e., very close to a linear dependence. Note that 
iRusseU et al] l|2013h ignore altogether the fact that the jet blob transitions from opti- 
cally thick to thin at light curve maximum, thereby missing a key piece of information. 
As a result, they do not have an analytic estimate of R and need to estimate tjct from 
the poorly constrained rise time to maximum of the radio flux. 



8 Conclusions and Future Prospects 

The continuum-fitting method has a number of virtues. A principal one is the simple 
and elegant model upon which it is based, namely the model of a thin, viscous accre- 
tion disk. This model was anti cipated a nd developed for black holes well before their 
existence was established fe.g.. iLvnden-Bell. 1969, 1 . Shortly thereafter, the analytic the- 
ory of thin disks was fully devel oped, an effort that culminated in the workhorse NT 
model l|Novikov fc Thorne|[l97i '). The most important predictions of this model have 



been validated recently via GRMHD simulations (Sections l4.2l and l5.2p . This venerable 
model, with the addition of a model of the disk's atmosphere (Sections 13.11 and 15. 3|) . 
measurements of three key parameters {D, i and M) and suitable X-ray data, allows 
one to estimate the inner-disk radius. Meanwhile, an abundance of strong observational 
and theoretical evidence allows one to identify this radius with the radius of the ISCO, 
which is simply related to the spin of the black hole. 

Another key virtue of the continuum-fitting method is an abundance of data for 
which the thermal disk component is strongly dominant. For most stellar black holes, 
one has many suitable archival spectra for which only a few percent or less of the 
thermal seed photons are upscattered by a hot corona into a Compton tail component 
of emission. In short, the continuum-fitting model is tried and true, and there is an 
abundance of available data for many stellar black holes obtained for the simplest and 
best understood state of an accreting black hole, namely, an optically-thick thermal 
disk. 

The spins and masses of ten stellar black holes are given in Table 1. Their spins span 
the full range of prograde values, and their masses range from 6 — 16 Mq. Setting aside 
the extreme spin of the transient GRS 1915-1-105, the three persistent black holes have 
higher spins and larger masses than their transient cousins. Furthermore, the high spins 
of these persistent, young black holes are unlikely to have been achieved via accretion 
torques, which implies that their spins are natal. The spins and masses of the ten black 



The ejection apparently shuts off, coincident ally, as the jet becomes optically thin. 
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holes in Table 1 provide their complete description and are the essential data for testing 
astrophysical models of how an accreting black hole interacts with its environment. 
They are likewise essential data for underpinning the physical theory of black holes, 
and for ultimately attempting tests of the no-hair theorem by, e.g., observing deviations 
from the multipoles predict ed by the Kerr metric, all of which are functions of just the 
two paramet ers a* and M l|Vigeland fc Hughesll201ol : Ijohannsen fc Psalti3l2O10l . l2Olll : 
iBambilboid ). 

To date, the most important application of the data in Table 1 is to one class of 
black hole jets, namely, those produced by a transient black hole near its Eddington 
limit as its spectrum switches from hard to soft. The peak radio luminosities of five of 
these transient microquasars correlate strongly with their spins, increasing by a factor 
of ~ 1000 as spin increases from ~ to > 0.95. Meanwhile, a simp le synchrotron 
j et mo del shows radio luminosity to be a good proxy for jet power. As lNaravan et al] 
(|2013l ) discuss in detail, the fitted relationship between jet power and spin (Figure [SJa) 



is not only a validation of the classic model of lBlandford fc Znaiekl (1973), it was also 
anticipated by GRMHD simulations showing that a jet can extract energy directly 
from a spinning black hole. 

During the next several years, one can hope to double the number of black holes 
with spins measured via the continuum-fitting method. It will be equally important 
to improve the quality of each measurement, largely by obtaining more accurate mea- 
surements of the parameters D, i and M, but also by making methodological advances 
and by pursuing more advanced GRMHD and atmosphere models of thin disks. The 
payoff for this effort will be the widening applications of these spin data to problems 
in astrophysics and physics. 

Especially important will be the possibility of validating both the continuum-fitting 
and Fe-line methods by comparing spin results obtained for individual stellar black 
holes. The validation of the Fe-line method is particularly important because it is the 
most direct approach to measuring the spins of AGN. Pr esently, concordant results 
are being obtained using these tw o leading methods (e.g., ISteiner et al.|[201ll . l2012"bl : 
iFabian et aLll2012l : lRevnoldsll2013l ). 

Two other methods for me asuring the spins of stellar black holes appear promising, 
namely, via X-ray polarimetry i Dovciak et al.|[2008l :lLi_et_al .l2009l:ISchnittman fc KrolikI 



^009) and high-freq uency quasiperiodic oscillations I T5r5^et''aLll2005 : RerniUarTTT^cClintoclj 



2006, : Belloni et al.l 12012, 1. The former method is stymied by a lack of data, and the 
latter by the lack of an adequate model. It is reasonable to hope that the mass and 
spin data in Table 1 will assist in identifying the appropriate physical model for the 
high-frequency QPOs. Because spin is such a critical parameter, it is important to 
attempt to measure it by as many methods as possible, as this will arguably provide 
our best check on the results. Stellar black holes are central to this effort because all 
of the methods of measuring spin mentioned above can be applied to them. 
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